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^ : Abstract 



(—] ' We present the fundamentals of geometrothermodynamics, an approach to study the properties 
of thermodynamic systems in terms of differential geometric concepts. It is based, on the one hand, 

QJ • upon the well-known contact structure of the thermodynamic phase space and, on the other hand, 

^ [ on the metric structure of the space of thermodynamic equilibrium states. In order to make these 

. ^ ! two structures compatible we introduce a Legendre invariant set of metrics in the phase space, and 

^^' demand that their pullback generates metrics on the space of equilibrium states. We show that 

, ^~ <. Weinhold's metric, which was introduced ad hoc, is not contained within this invariant set. We 

psj ' propose alternative metrics which allow us to redefine the concept of thermodynamic length in an 

>■ 
Tij" ' invariant manner and to study phase transitions in terms of curvature singularities. 

■^ ■ PACS numbers: 05.70.-a, 02.40.Ky 
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I. INTRODUCTION 

One of the first applications of differential geometry in thermodynamics is due to Gibbs^ 
and to Caratheodory^. Those studies were further developed in the seventies specially by 
Hermann:^ and later by Mrugalai^i^, leading to a geometric approach in which the thermody- 
namic phase space and its natural contact structure are the basic elements of the construc- 
tion. Weinhold^ proposed a second approach with a metric defined ad hoc as the Hessian 
of the internal energy. This approach has been intensively used to study the properties of 
the space generated by Weinhold's metrio^i^, the thermodynamic lengtb^ii^iiii, the physical 
properties of various two-dimensional thermodynamic systemsi^'i^ii^ii^'i^, and the associated 
Riemannian structure^»i2ii^. In an attempt to understand the concept of thermodynamic 
length, Ruppeiner— introduced a metric which is conformally equivalent to Weinhold's met- 
ric. The multiple applications of Ruppeiner's metric have been reviewed in=i and more 
recent results are included ii>iii22ii^i2i. These two approaches have been the subject of sep- 
arate analysis since the seventies, but the question about their relationship remains open. 
Several attempts were carried ouir^^, resulting in a metric which was non-invariant under 
Legendre transformations. This unpleasant property means that a given thermodynamic 
system has different properties, when different thermodynamic potentials are used. 

The main purpose of the present work is to show that both approaches can be unified 
into a single approach by using purely mathematical considerations. We will call this ap- 
proach geometrothermodynamics. The idea is simple. On the phase space we first consider 
an arbitrary metric which, by means of the pullback, induces a metric on the space of equi- 
librium states. Then we derive the conditions for these metrics to be invariant with respect 
to arbitrary Legendre transformations, making the contact structure compatible with the 
Riemannian structure. Geometrothermodynamics is a unifying approach which allows us to 
treat thermodynamics in a geometric language either at the level of the phase space or in the 
space of equilibrium states. The first attempt towards a formulation of geometrothermody- 
namics was presented in^ where the connection, torsion and curvature of the fundamental 
Gibbs 1-form were used to construct a geometry for equilibrium states. 

This paper is organized as follows. In Section [TTl we present the main geometric struc- 
tures which are necessary for the formulation of geometrothermodynamics. This includes 
the contact structure, Legendre invariance, and the Riemannian structures. In Section IIIII 



we present the special case of two-dimensional thermodynamic systems. We analyze Wein- 
hold's approach and show that it does not satisfy all the requirements necessary to describe 
thermodynamic systems in an invariant manner. We find the conditions for a metric to be 
Legendre invariant and show the existence of metrics with this property. Furthermore, in 
Section HVl we analyze the simple case of an ideal gas and propose a metric whose curvature 
can be used as a measure of thermodynamic interaction. The same metric is used to pro- 
pose an invariant definition of thermodynamic length. In Section |V] we analyze the van der 
Waals gas and prove that critical points and phase transitions of any two-dimensional ther- 
modynamic system are characterized by curvature singularities of a specific metric. Finally, 
Section |VT] is devoted to discussions of our results and suggestions for further research. 

II. FUNDAMENTALS OF GEOMETROTHERMODYNAMICS 

A. Contact structure 

The first step to introduce the language of differential geometry in thermodynamics is the 
definition of the (2?2-|-l)-dimensional thermodynamic phase space T. It can be coordinatized 
by the set Z^ = {$,£'",/"}, where $ is the thermodynamic potential, and E'^ and 1°" 
denote the extensive and intensive variables, respectively. Here we adopt the conventions 
A = 0, ...,2n and a = 1, ...,n so that $ = Z°, E"- = Z"-, and /'^ = Z"+". Furthermore, we 
introduce the fundamental Gibbs 1-form (summation over repeated indices) 

e^ = rf$ - d^J^dE' , 6,t = diag(l, 1, ..., 1) (1) 

The pair (T, O^) is called a contact manifold^ if T is differentiable and 0^ satisfies the 
condition B^ A (rfO^)" ^ 0. Consider now the n— dimensional space S spanned by the 
coordinates E"-. This can be realized by means of the smooth mapping ip : S — > T 

^■.{E'')^^{<^,E-,n , (2) 

with $ = ^{E"-). We define the space of thermodynamic equilibrium states S as the subspace 
of T given by the embedding mapping (p : S ^ T for which the condition 

ip*{e^)=ip*{d^-5abPdE')=0 (3) 



holds, where ip* represents the pullback. This imphes the relationships 



'* «./' (4) 



which are the standard conditions of thermodjTianiic equilibrium. Consequently, on the 
space of equilibrium states £^ condition ([3]) leads to the first law of thermodynamics 

d^ - Sabl^dE^ = . (5) 

To distinguish one thermodynamic system from another in the space of equilibrium states, 
one can specify the fundamental equation^l which, in the representation used here, is con- 
tained in the embedding mapping ([2]) through the relationship $ = $(i5"). This construction 
must be complemented with the second law of thermodynamics 

dE^^^ - ° ' (^) 

which is also known as the convexity conditioii^^i^S. The thermodynamic potential must 
satisfy the homogeneity condition ^{XE"") = X^^^E"") for constants parameters A and (3. 
Differentiating this homogeneity condition with respect to A and evaluating the result at 
A = 1, we get an expression which transforms into Euler's identity 

P<!>{E-) = 6aj'E- , (7) 

after using the duality condition (jl]). Calculating the exterior derivative of Euler's identity 
and using ([5]), we obtain the generalized Gibbs-Duhem relation 

(1 - P)SabrdE' + SabE'^dl' = . (8) 

The classical expressions for Euler's identity and Gibbs-Duhem relation are obtained from 
the above equations by putting /3 = 1. 

B. Legendre invar iance 

Legendre transformations are a special case of contact transformations which leave in- 
variant the contact structure of T. In physical terms, Legendre invariance means that the 
thermodynamic properties of a system are independent of the thermodynamic potential used 
to describe it. Let us consider a partial Legendre transformation^ 

{Z^} — > {Z^} = {$, E", r} (9) 



$ = $ - SkiE'^P , E' = -r, E^ = E^, r = E\ P = P , (10) 

where i U j is any disjoint decomposition of the set of indices {1, ..., n}, and k,l = 1, ..., i. 
In particular, for i = {1, ..., n} and i = we obtain the total Legendre transformation and 
the identity, respectively. In the further discussions we will always consider total Legendre 
transformations for the sake of simplicity. The main results will not be affected by this 
simplification. 

The invariance of the contact structure introduced above becomes evident from the fact 
that in the new coordinates {Z^} the Gibbs 1-form ([1]) becomes 6^ = d^ — 5abI°'dE^. A 
straightforward calculation proves also the invariance of the duality condition (jlj) and the 
first ([5]) and second (Ej) laws of thermodynamics. Notice that a Legendre transformation 
interchanges the character of the extensive and intensive variables. For this reason, to 
establish the validity of Euler's identity and, consequently, the Gibbs-Duhem relation in the 
new variables it is necessary to identify first the extensive variables E'^ and the fundamental 
equation $ = $(i5"). Then, Legendre invariance follows from the fact that the duality 
condition (j4]) and the first law ([5]) are invariant. 

C. Metric structure 

Let us consider a non-degenerate metric G which induces a Riemannian structure on the 
thermodynamic phase space T. Now we can formulate the main statement of geometrother- 
modynamics. A thermodynamic system is described by a metric G which is called a ther- 
modynamic metric if it satisfies the following conditions: 

1. G is invariant with respect to transformations which do not modify the contact 
structure of T . In particular, G must be invariant with respect to Legendre transformations. 

2. G induces in the space of equilibrium states £ an invariant metric g by means of the 
mapping 

^\G)=g. (11) 

The idea behind condition 1 is that a thermodynamic metric must describe a thermodynamic 
system, independently of the coordinates used in T. This condition is necessary in order 
for geometrothermodynamics to be able to describe thermodynamic properties in terms of 
geometric concepts in a manner which must be invariant with respect to changes of the 
thermodynamic potential. Condition 2 establishes the relationship between the geometry 



on T and the geometry on £ by using the same tools that are used to define equihbrium 
states in the contact structure approach. 

It is clear that the Legendre invariance of g follows from the invariance of G, if they are 
related by Eq. fllll) . In fact, if we denote by (5 = Gi^Z"^) the metric obtained from G by 
applying the Legendre transformation fllUp . and define G'{Z^) = G{Z^ = Z^), then the 
Legendre invariance of G can be reached by demanding that G{Z^) = G'{Z^). Furthermore, 
since g = (p*{G) and g' = (p*{G') with (p defined as in Eq.([2]) with Z^ = Z^, then we have 
that g = 'f*{G) = (p*{G') = g', i.e. g is an invariant metric. 

Under a total Legendre transformation, Z^ — > Z^, as given in Eq. lTTOl) . the components 
of the metric G transform as 



G 



AB 



~ _dZ^dZ^ 

\jr AR — — ~ = — ^CD 

dz^dzB 



(12) 



with the transformation matrix 
/ 1 -Z"+i 



dz^ 














1 



















1 



-z^ 
-1 








-Z"\ 










-1 





(13) 



yO ... 1 ... / 

where the indices A and B represent rows and columns, respectively. Since the determinant 
of this matrix is equal to one, the inverse transformation exists, and it determines the 
diffeomorphism corresponding to a total Legendre transformation. 

To guarantee that the metric G in T is Legendre invariant it is sufficient to demand that 
each of its components preserves its functional dependence after a Legendre transformation. 
Let Gab{,Z^) = Gab{,Z'~^ -^ Z^) represent the components of the metric G obtained after 
a Legendre transformation (TTOl) of the metric G, and let G^^(Z'") = Gab{.Z'~' = ^'") be the 
components of the metric G' obtained by replacing each coordinate Z^ by its counterpart 
Z^ in the metric G (no Legendre transformation is applied). Then, if we demand that the 
condition 

Gab{Z^) = G'abCz"") (14) 
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holds for each value of A and B, we guarantee that the metrics G on T as well as g on S 
are Legendre invariant. 

An important issue to be addressed when imposing invariance conditions on metrics is 
the one related to the existence of solutions. In our case, the point is whether there exist 
metrics which satisfy the condition of invariance under Legendre transformations as given, 
for instance, in Eg. (1141) . In fact, it is easy to construct an invariant (2n + 1)— dimensional 
metric by defining the Gibbs metric as the "square" of the Gibbs 1-form ([1]), i.e. 

G = e^ ® e^ = e^ = rf$2 - 25abI'"dE^d^ + SabScdl'^I^dE^dE'^ . (15) 

The invariance of this metric follows from the invariance of the Gibbs 1-form. Since the Gibbs 
metric does not possess components proportional to t/J", its determinant vanishes. How- 
ever, it contains an {n + l)-dimensional metric in the sector {(i$, dE*^} with non-vanishing 
determinant. Despite its degeneracy, later on we will use the Gibbs metric metric to con- 
struct explicit examples of non-degenerate 5-dimensional metrics which satisfy the Legendre 
invariance condition. 

III. SYSTEMS WITH TWO DEGREES OF FREEDOM 

A. Two-dimensional geometrothermodynamics 

Let us consider a mono- component thermodynamic system with only two degrees of 
freedom. The thermodynamic phase space T is therefore a 5-dimensional manifold for 
which we choose the coordinates 

{Z^} = {$, E^ n = {U, S, V, T, -P} , (16) 

in accordance with standard notations of thermodynamics^^. The Gibbs 1-form 

e^ = dU- TdS + PdV (17) 

defines the contact structure of T. The 2-dimensional space of equilibrium states S is 
determined through the embedding mapping 

y. : {S, V) ^ [U{S, V), S, V, T{S, V), P{S, V)] , (18) 



so that (f*{Q^) = implies the first law of therinodynainics dU = TdS — PdV with 

-P . (19) 



dU _ dU 

'dS~'dV 



Notice that the fundamental equation U = U{S, V) is contained in the definition of the 
mapping flTSl) . as mentioned in Section [II A[ The dependence of the intensive variables T 
and P in terms of the extensive variables S and V is induced through the conditions (fT9l) . 
The Euler identity and Gibbs-Duhem equation in this case can be written as 



PU-TS + PV = , SdT - VdP + (1 - (3)dU = 



(20) 



and reduce to the usual form when /? = 1, i.e., when 5" and V are strictly extensive variables. 
A total Legendre transformation {U, S, V, T, —P) i — > {U , S, V, T, —P) is defined in this 
case by means of the relationships 



U = U-ST + PV 



S 



-T,V = P,T = S,P 



-V 



(21) 



and the corresponding transformation matrix reduces to 

/ 1 -f P -S v\ 



dZB 












-1 

















1 





1 

















-1 





) 



(22) 



When compared to the general transformation matrix flT3l) . this matrix presents several 
changes in the signs which are due to the choice of the pressure as P = —Z^. Also, it is 
convenient to write the general 5-dimensional metric of T as {I, J, ... = 1, ...4) 



G = GoodU'^ + 2GojdUdX^ + GijdX^dX^ , X^ = {S, V, T, P) 



(23) 



so that the negative sign of Z^ does not affect the sign of the metric components. 

According to our construction of geometrothermodynamics, we must demand that the 
metric fl25]) be invariant with respect to Legendre transformations. If the components Gab 
are given as functions of the coordinates of the thermodynamic phase space, a total Legendre 
transformation leads to a set of algebraic equations which are explicitly given in Appendix 
lAt where we also show the existence of solutions. For a general 5-dimensional metric G we 
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can also compute the corresponding 2-diniensional metric g by using the pullback of the 
mapping (fTSj) . i.e. 

g = ^*{G) = gndS^ + 2gudSdV + g22dV^ , (24) 



/2 



(25) 



9ii = G*!! + 26*01^ + GqqT + Gs^Ugg + G/iJJ^ 

+ 2[/55'(Gi3 + GoS^ ^ G'34f/sv') ■" 2f/5'y(G'i4 + 6*04^) , 

5'12 = Gi2 — P{Gqi + GqqT) + G'o27' — ^34^/5^ 

+ ^5^(^13 + Gos?" — G24 + G04-P + G33f/5S + Gi/JJvv) 

+ Uss{G23 — G03-P — G34f/yy) — f/yy(Gi4 + ^047") , (26) 

922 = G22 — 2G02-P + GoqP + GssUgy + GiiUyy 

+ 2Uvv{-G24 + G04P - G34f/5V) + 2Usv{G23 - G03P) , (27) 

where the subscript of U represent partial derivatives, and we used equations flT9|) in the 
form Us = T, Uy = — -P- This is the most general metric that can exist on the space of 
equilibrium states £. The further specification of the fundamental equation U = U{S, V) 
fixes completely the structure of g in terms of the components Gab which, in general, 
become arbitrary functions of 5* and V, due to the definition of the embedding mapping 
(ITSl) . This arbitrariness can be restricted by considering only Legendre invariant metrics 
G. In fact, if we apply a total Legendre transformation on the 5-dimensional metric (l23l) . 
the new metric G = G{Z^ -^ Z^) will be Legendre invariant if its components satisfy the 
conditions (]A3p . Furthermore, if we start from a Legendre invariant metric G and perform 
the pullback determined by means of the mapping 



ip{S, V) ^ [U{S, V), S, V, T{S, V), P{S, V)] , (28) 

we will obtain a new metric g = (p*{G) whose components satisfy the relationship 

-gabCS. V) = gUS, V) = ga,{Z^ = ^^, Gab = Gab) , (29) 



where gab are given in Eqs. fl25|) - fl27j) . In other words, if the starting metric G is Legendre 
invariant, the resulting metric g = f*{G) is also Legendre invariant and its components are 
given again as in Eqs .(1251) -( 1271) . 



We found the explicit form of the most general 2-dimensional metric g on S, starting 
from a general 5-dimensional metric on T. The set of algebraic equations fl25|) - fl27|) also 
allows us to analyze the inverse problem, i.e. if we know the explicit components Qab of a 
metric on S, we can find from equations fl2^ - fl7r|) the most general form of the metric G 
which generates g by means of the puUback. In this manner, we can determine whether a 
given metric g can be obtained from a Legendre invariant metric G. 

As a simple example let us consider the 2-dimensional Euclidean metric 

g^ = dS^ + dV^ . (30) 

Introducing these values into fl25|) - fl27|) . we see that the generating metric can be written as 

G^ = GQo[dU^ + {TdS - PdVf] + dS^ + dV^ + 2GQidX^Q^ , (31) 

where X^ = {S, V, T, P) and 9^ = dU — TdS + PdV is the fundamental Gibbs 1-form for 
systems with two degrees of freedom. The components Gqa remain as arbitrary functions of 
the coordinates of T. This is an example of the freedom which exists in the determination 
of 2-dimensional metrics on S in terms of the pullback of 5-dimensional metrics on T. A 
straightforward calculation shows that the determinant of G^ vanishes identically, making 
this metric unsuitable for describing a Riemannian structure on T. Nevertheless, to show 
the details of the procedure, we want to explore the Legendre invariance of the Euclidean 
metric on S. Applying a total Legendre transformation (|2T1) to G^ , we obtain 

G^ = GqqQI + df^ + dP^ 

+ 2{-Goidf + Go2dP + GosdS - GoAdV)e^ , (32) 

where we used the Gibbs-Duhem relation VdP — SdT = in the new variables to simplify 
the final result. Finally, we calculate the pullback determined by the embedding mapping 
(1281) and obtain 



~g^ = <^*(G^) = ( S^] dE^dE' 

i^h + Uly)dS' + 2U^y{Uss + Uyy)dSdV + {illy + Uly)dV^ , (33) 
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an expression which is functionally different from the starting metric (130|) . This shows that 
the 2-dimensional Euclidean metric g^ is not invariant with respect to Legendre transfor- 
mations and, consequently, is not appropriate for describing any thermodynamic system in 
the context of geometrothermodynamics. 

B. Weinhold's approach 

Weinhold proposed to study the geometry of thermodynamic systems by means of a metric 
which is introduced ad hoc on the space of equilibrium states £. In the case of a system with 
two degrees of freedom, it corresponds to a 2-dimensional metric with components given in 
terms of the second derivatives of the internal energy with respect to the extensive variables. 
When the entropy representation is used, an alternative metric can be defined again as the 
Hessian of the entropy (Ruppeiner's metric) which, however, turns out to be conformally 
equivalent to Weinhold's metric. Since any metric on a 2-dimensional space is conformally 
equivalent to the Euclidean metric, it seems reasonable to expect a conformal equivalence 
between all 2-dimensional metrics, when arbitrary diffeomorphisms are allowed. 

In the energy representation we are using here, Weinhold's metric is given by 

g"^ = (g§^) dE'^dE" = UssdS^ + 2UsvdSdV + UyvdV^ . (34) 

An important property of Weinhold's metric is that it is positive definite since its compo- 
nents coincide with the matrix components of the second law of thermodynamics ([6]). To 
determine the 5-dimensional metric on T that generates Weinhold's metric, we introduce 
its components into equations fl2^ - fl27|) . After an algebraic rearrangement, the resulting 
metric can be written as 

G^ = Goo[dU^ - {TdS - PdV)^] + dSdV - dTdP + 2GojdX^Q^. (35) 

We see that the generating metric is not unique, but it possesses five degrees of freedom which 
are contained in the five arbitrary functions Gqa- The rest of the components of Gab are 
fixed through the components of Weinhold's metric. A special case of this generating metric 
was proposed in-^^ and corresponds to the particular choice Goo = l? Coi = ~T, Gq2 = 

-Pj G^03 = Gqa, = 0. 

The next important question is whether the metric G^ is invariant with respect to 
arbitrary Legendre transformations. Introducing the explicit components of G'^ into the 
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conditions of Legendre invariance given in Eqs. (lA3p . one can see that they are not satisfied, 



independently of the values of the arbitrary functions Gqa- Consequently, Weinhold's metric 
is not Legendre invariant, making it inappropriate for describing thermodynamic systems, 
because the properties of the system would depend on the thermodynamic potential used 
in g^ . This, of course, is a disadvantage which is sufficient for dismissing this metric 
as unphysical. Nevertheless, to find the cause of the problem, we perform a total Legendre 
transformation on G^ to obtain a new metric G^ . Then we calculate its puUback ip*{G'^) = 
g^ and obtain the metric 

~gW = _ (^E\ dE-d& , (36) 

\dE^dE^J ^ ' 

which obviously contradicts the second law of thermodynamics (|6]). We see that the non- 
invariant character of this metric lies in the change of the metric signature. One could 
naively think that the additional "parity transformation" U —>■ —U would preserve invari- 
ance. However, this parity transformation on S is not allowed because it does not correspond 
to a Legendre transformation on T. Another possibility could be to consider a "quadratic" 
Weinhold's metric so that the signature becomes invariant. Unfortunately, this simple al- 
ternative is also not allowed because the resulting metric is not invariant. In fact, this is 
the example we presented in Section IIII A[ There, it was shown that the Euclidean metric 
g^ turns into the quadratic Weinhold's metric ( l33l) after a Legendre transformation. The 
inverse transformation would then transform the quadratic Weinhold's metric into the Eu- 
clidean one which, as we mentioned above, is not suitable for a geometrothermodynamic 
approach. 

C. Legendre invariant metrics 

Following the procedure described in Appendix [Al it is in principle very simple to find 
Legendre invariant metrics. By assuming a polynomial dependence of the metric components 
in terms of the thermodynamic variables, we derive in Eqs. (]A4|) a particular solution which 
possesses twelve arbitrary real constants. Clearly, more general solutions can be obtained 
by starting from metric components with a different functional dependence. 

In this Section, however, we will derive several Legendre invariant metrics by using a 
simple and intuitive method. From the definition of a Legendre transformation given in fITOl) . 
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it follows that metrics with no 0?$ terms and combinations of the form YlAi^-^'^)^ + (c/J")^] 
and Ylai^'^^'^'^ + E^-dl"-]"^ are invariant with respect to Legendre transformations. This can 
be used as guide to construct invariant metrics. In particular, for systems with two degrees 
of freedom let us consider the following metric 

G^ = 92 + ciTSdTdS + CiPVdPdV + dS'^ + dV^ + dT"^ + rfP^ (37) 

which is Legendre invariant for any values of the real constants Ci and C2. The term 0^ has 
been added in order to avoid the degeneracy of G^ . The pullback of this metric ip*{G^) = g^ 
generates a 2-dimensional metric on S which in a matrix representation can be written as 

yiiciTS ~ C2PV)Usv -C2PVUVV J 

where 1 is the unit matrix and g^ represents Weinhold's metric flMl) . The special case 
Ci = C2 = corresponds to Hernandez-Lacomba metric^^ which is also Legendre invariant. 
Consider now the following 5-dimensional non-degenerate metric 

Qii = 92 + {TS - PV){dSdT - dVdP) , (39) 

which is also Legendre invariant. A straightforward computation of its pullback yields 

g^^ = U {UssdS^ + 2UsvdSdV + UyvdV^) (40) 

where we used Euler's identity (l20l) . We will see in the next Section that this metric resembles 
the behavior of Weinhold's metric in the case of an ideal gas. 

Clearly, the metrics presented above can easily be generalized to the case of a system 
with n degrees of freedom, preserving Legendre invariance. Dropping constant factors, the 
generalized (2n + l)-dimensional metrics can be written as 

G^ = Ql + SabE'TdE^'dl'' + 5ab{dE''dE^ + dl^dl^) , (41) 

C' = Ql + {dabE'^l') {d.dE'dl") , (42) 

where 0^ is the Gibbs 1-form ([T]) for a thermodynamic system with n degrees of freedom. 
All the examples we have presented here are characterized by the fact that the metric 
components can be written as polynomials of the thermodynamic variables. Solving the 
algebraic conditions described in Appendix |X1 it is also possible to derive Legendre invariant 
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metrics of T whose components are rational functions of the variables. The corresponding 
pullback generates metrics on £ with rather cumbersome structures which, nevertheless, are 
allowed in the context of geometrothermodynamics. 

IV. THE IDEAL GAS AND THERMODYNAMIC LENGTH 



In our opinion, one of the most interesting properties of Weinhold's metric is that its 
curvature vanishes in the case of a non-interacting thermodynamic system, i.e. an ideal 
gas. For more general systems, like the van der Waals gas, the curvature turns out to be 
non-vanishing and, therefore, it has been postulated that it could be used as a measure 
for thermodynamic interaction. Although we have shown that Weinhold's metric is not 
compatible with the requirements of geometrothermodynamics, we would like to preserve in 
our approach this interesting interpretation of curvature. To this end, we need to show that 
there exists at least one Legendre invariant metric whose curvature vanishes in the case of 
an ideal gas. We will show that in fact we already have a metric with this property. 

Let us recall that for a 2-dimensional metric gab there is only one independent component 
of the Riemann curvature tensor Rabcd which completely determines the curvature scalar 
R = g"''^g^'^ Rabcd- in turn, the curvature scalar can be expressed as follows 



R 



1 



Vdet{g) 



9n,2 — 9i2, i 



922^1^2.912^ 

Vdet((7) 



2detigy 



det{H) 



(43) 



with 



/ 



H 



\ 



(44) 



9ll 9l2 922 

9ll,l 9l2,l 922,1 

\ 9ll,2 9l2,2 922,2 I 

where a comma denotes partial differentiation. As mentioned before, all the information 
about a specific system is contained in the fundamental equation. In the case of an ideal 
gas it can be written as2i 

U{S,V) = V-'/'exp(^^^ , (45) 

where /c is a constant. From this equation one can calculate the standard thermodynamic 
variables which are related through the equation of state PV = kT. Finally, the constant 
entropy is also determined by means of Euler's identity. Introducing this fundamental equa- 
tion into the components of the metric g^ given in (|38|) and computing the corresponding 
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curvature scalar, we obtain a rather cumbersome expression which is not a constant and 
is always different from zero, independently of the values of ci and C2. We conclude that 
the curvature of the metric g^ is non zero for an ideal gas and, consequently, it cannot be 
used as a measure of thermodynamic interaction. Consider now the metric g^^ as given in 
Eg. (1401) . A straightforward computation of its components with the fundamental equation 
(H5|) leads to the following result: 

gl'g = ^^^{2V^dS^ - AkVdSdV + hk''dV^) . (46) 

It is then easy to show that the corresponding scalar curvature vanishes identically for 
an ideal gas. Moreover, for more general systems, like the van der Waals gas, the scalar 
curvature of g^^ does not vanish (see below). This is exactly the property we need in order 
to introduce the concept of thermodynamic length in an invariant manner. 

Let us consider a system with n degrees of freedom whose phase space T is described 
by the general metric G^^ as given in ( l42l) . The computation of its puUback, using the 
conditions of thermodynamic equilibrium (j4j) and Euler's identity (JTj), yields 

/ 52$ \ 

^^^ = $ -— -— dE'^dE^ . (47) 

This metric is defined on the space of equilibrium states £ and, in contrast to Weinhold's 
metric ( IMj) . is Legendre invariant due to the presence of the conformal factor $. Since the 
curvature of this metric is a measure of thermodynamic interaction, it seems plausible to 
use it also as a measure of "length" on £^. If ti and ^2 represent two different states on £, 
we define their thermodynamic length as 



i2 



(92$ 
$ ( — ) dE^dE' 
dE^dE^ 



1/2 



(48) 



Similar definitions have been proposed by Weinhold^ and Ruppeiner— , using, however, non 
Legendre invariant metrics. We believe that our proposal fHHj) will be useful for clarifying 
the true significance of thermodynamic length. 

V. THE VAN DER WAALS GAS AND PHASE TRANSITIONS 

The fundamental equation for the van der Waals gas can be written as 

2S 
P 3k CI 
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where a, b, and k are constants. The thermodynamic interaction is determined through the 
constants a and b. For this reahstic gas we find that the metric (HOj) reduces to 



9",w = ^U{U + a/V) 



2 



2dS' - J^dSdV + jJ^dV' 
V — b [V — b)^ 



2aU 
17^ 



dV' (50) 



and the corresponding scalar curvature can be expressed as 



aV{U,V,a,b) 
U^{PV^-aV + 2aby ^ ' 

where P(f/, V, a, 6) is a polynomial which is always different from zero for any real values of a 
and b. In the limiting case a = 6 = 0, the curvature vanishes and we turn back to the case of 
an ideal gas. The curvature also vanishes for a = and & 7^ 0, indicating that the constant 
a is responsible for the non-ideal thermodynamic interactions, whereas the constant b plays 
a qualitative role in the description of interactions^. When the condition 

PV"^ - a1^ + 2a6 = (52) 

is satisfied, there are true curvature singularities which could be interpreted as showing the 
limit of applicability of geometrothermodynamics. Indeed, condition fl52|) for the van der 
Waals gas corresponds to the limit of thermodynamic stability^ from which the information 
about critical points and phase transitions can be obtained. Notice that condition fl52|) for 
curvature singularities represents the failure of local stability and, consequently, contains 
information about second order phase transitions only. No information about first order 
phase transitions can be obtained in this case from the analysis of curvature. 

The above result can be formulated independently of the thermodynamic system as fol- 
lows. The expression for the scalar curvature (jH]) shows that it is well defined only if 
det((7) > and becomes singular if det(5') = 0. In the case of the metric g^^ given in 
Eq. fHUl) . these two conditions lead to 

VssUvv - t/|^ > , (53) 

which corresponds to the convexity and stability conditions for any thermodynamic system 
with two degrees of freedom in the energy representation^^. We conclude that the curvature 
singularities of the metric g^^ on E contain the information about the phase transitions of 
the thermodynamic systems. 
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VI. CONCLUSIONS 

In this work we formulated the fundamentals of geometrothermodynamics which unifies 
the contact structure of the thermodynamic phase space T with the metric structure of the 
space of thermodynamic equilibrium states £ d T . The unifying object is a Riemannian 
metric structure on T which must be invariant with respect to Legendre transformations, 
and compatible with the metric structure induced on £ by means of the pullback. It was 
shown how the first and second laws of thermodynamics, as well as Euler's identity and the 
Gibbs-Duhem relation, are invariant with respect to Legendre transformations. We found 
explicitly the conditions under which the Legendre transformations are also isometrics for 
a metric structure on T . This leads to a system of algebraic equations for the components 
of the metric. Several particular solutions were found for this system. All the solutions 
found turned out to be represented as polynomials in terms of the coordinates of T, but it 
is possible to derive more general solutions. 

We analyzed in detail the case of thermodynamic systems with two degrees of freedom. 
Weinhold's approach was investigated for this case and we found that it does not satisfy the 
conditions of Legendre invariance and, therefore, cannot be used to describe thermodynamic 
systems in an invariant way. In particular, we saw that it is in conflict with the second law 
of thermodynamics. As explicit examples of simple thermodynamic systems we investigated 
the ideal gas and van der Waals gas. Within the set of Legendre invariant metrics included 
in this work, we found one particular metric whose curvature vanishes in the case of an ideal 
gas, and is different from zero otherwise. This result makes this metric especially suitable 
for measuring the interaction in thermodynamic systems. 

We proposed a definition of thermodynamic length which, in contrast to other definitions 
known in the literature, is Legendre invariant. Furthermore, we proved that in the case of 
thermodynamic systems with two degrees of freedom, the curvature singularities of a specific 
Legendre invariant metric indicates the existence of critical points and phase transitions. 

It would be interesting to continue the study of the metrics presented here in the case 
of more complicated thermodynamic systems like the Ising model, multi-component ideal 
gas, etc. This would shed light on the significance of these metrics. On the other hand, our 
approach opens the possibility of investigating simultaneously the geometric structure not 
only in the space of equilibrium states, but also at the level of the thermodynamic phase 
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space. These tasks will be treated in future works. 

Our geometrotherniodynamic approach opens new possibilities to study thermodynamic 
systems in terms of geometric objects in an invariant manner. At the same time, we will have 
to face new challenges. In particular, we have found a huge arbitrariness in our description. 
For instance, the study of the ideal gas allowed us to stand out the metric g^^ given in 
Section HV] because its curvature vanishes in the limiting case of an ideal gas. However, our 
results seem to indicate that this is not a special property of that metric. With the help of 
the algebraic conditions presented in Appendix [Al it would be possible to find more metrics 
with the same property. It seems that we need additional criteria to select or classify metrics 
on the thermodynamic phase space. Understanding the metric in this context represents an 
interesting problem for future investigations. 
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APPENDIX A: LEGENDRE TRANSFORMATIONS FOR SYSTEMS WITH 
TWO DEGREES OF FREEDOM 

The components Gab of an arbitrary metric in a 5-dimensional thermodynamic phase 
space transform under a Legendre transformation flTOj) according to the rule 



Gr/iR — S> CjrAn — — = — Gen 1 (Al) 

where the transformation matrix is given in Eq. fl22l) . Using the representation fl23|) . a total 
Legendre transformation leads to 



18 



G^oo = G^oo 5 

Goi = Go3 — GoqT , 

Go2 = —Go4 + GoqP , 

Go3 = ~{Goi + Goo 5') ) 

Gq4 = Go2 + GooV^ , 

Gil = G^33 ~ SGosT + Goo^ > 

Gl2 = —Gs4 + G04T + G03P — GqqPT , 

Gi3 = — Gi3 — GosiS* + GqiT + GqqST , (A2) 

Gi4 = G23 + G03V — G02T — GoqTV , 

G22 = G44 — 2G04P + GqoP , 

G23 = Gi4 + Go4'S' — GqiP — GqqPS , 

G24 = —G24 — G()4V + G02-P + Goo-PV^ , 

G33 = Gil + 2Goi'S' + Goods' , 

G34 = — (G12 + Go25' + GoiV^ + Goo'S'K) , 

G44 = G22 + 2002^ + GooV' • 

This represents a system of algebraic equations for the components Gab which, if fulfilled, 
guarantees Legendre invariance. The explicit dependence of the above system from the 
coordinates of the thermodynamic phase space points to a particular solution which can 
be derived by considering an ansatz where Gqa depends linearly on Z"^ and the GikS are 
quadratic functions of the coordinates. A lengthly but straightforward computation shows 
that the metric 
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Goo 


= Ci , 


Goi 


= C2T, 


Go2 


= -C2P, 


Go3 


= -{Ci + C2)S , 


GoA 


= (Cl+C2)y, 


Gil 


= CsT^ + C4 , 


Gl2 


= -C5TP + ce , 


Gi3 


= C7ST , 


Gi4 


= -CsVT + Cg , 


G22 


= CioP' + Cii , 


G23 


= -CsPS - eg , 


G2A 


= C,2VP , 


G33 


= (ci + 2c2 + 03)5^ + C4 , 


G34 


= -{c, + 2c2 + C5)SV-c, 


G44 


= {Ci + 2C2 + Cio)V^ + Cu 



(A3) 



is completely invariant with respect to total Legendre transformations. Here Ci, ...,Ci2 rep- 
resent arbitrary real constants and in general det{GAB) 7^ 0. This particular solution shows 
that there exist non-trivial thermodynamic metrics which can be used to introduce Rieman- 
nian structures in the thermodynamic phase space. It is possible to derive other particular 
solutions by assuming a polynomial dependence of higher degree for all metric components 
Gab- We also investigated the case of rational and exponential functions, and obtained 
similar results. 
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